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ABSTRACT 

We describe the number counts and spatial distribution of 239 Distant Red Galaxies (DRGs), 
selected from the Early Data Release of the UKIDSS Ultra Deep Survey. The DRGs are iden- 
tified by their very red infrared colours with (J — K)ab > 1-3, selected over 0.62 deg'^ to 
a 90% completeness limit of Kab — 20.7. This is the first time a large sample of bright 
DRGs has been studied within a contiguous area, and we provide the first measurements of 
their number counts and clustering. The population shows strong angular clustering, inter- 
mediate between those of A'-selected field galaxies and optical/infrared-selected Extremely 
Red Galaxies. Adopting the redshift distributions determined from other recent studies, we 
infer a high correlation length of tq 11 h^^Mpc. Such strong clustering could imply that 
our galaxies are hosted by very massive dark matter halos, consistent with the progenitors of 
present-day L > L^, elliptical galaxies. 
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1 INTRODUCTION 

A new near-infrared selection technique has been developed in re- 
cent years to sample galaxies in the high-redshift Universe. By re- 
lying on purely near-infrared colours, this potentially avoids many 
biases which are inherent in optical te chniques, particula rly for 
detected dusty and/or evolved galaxies. iFranx et al] j2003h argue 
that the simple (J — K)ab > 1-3 colour selection criteria pro- 
duces a sample that is mainly populated by galaxies at 2: > 2, 
at least at faint A"-band magnitudes (Kab 21). These are the 
so-called Distant Red Galaxies (herea fter DRGs). In the Faint In- 
frared Extragalactic Survey (FIRES), iFranx et al] ( |2003|) selected 
14 candidate galaxies at z > 2 to a depth of Kab < 24.4, 
of wh i ch 6 were spectro scopically confirmed jvan Dokkum et all 
l2003h . lLabbe et alj j2005h found that approximately 70% of DRGs 
are dusty star forming galaxies and the remaining 30% are pas- 
sively evolved galaxies. 
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Work by iRudn ick et all j2003l) suggest that DRGs may be a 
signi ficant constituent of the z ~ 2 — 3 universe in terms of stellar 
mass. lForster Schreiber et all j2004h demonstrated that the average 
rest-frame optical colours of DRGs fall within the range covered by 
normal galaxies lo cally, unlike the Lyman-break galaxies (LBGs - 
ISteidel et aljT99^ which are typically much bluer Larger samples 
of DRGs covering a wide range in stellar mass are now required 
to fully understand the importance of this population. In particular, 
studies conducted so far have (by necessity) concentrated on DRGs 
selected over relatively small areas, and very little is known about 
the bright end of this population. 

In term s of st ellar mass, metalicity and star formation rate 
iReddv et al.l j2005h find strong similarities between optically- 
selected and near-infrared selected galaxy samples. Clustering of- 
fers an alternative way to to study these populations. At large 
scales, the galaxy distribution is dominated by dark mater halo 
clustering, which is a strong function of halo mass. Several studies 
have measured strong clustering strength for high reds hift galaxies 
selected in the near-infrared (ro = 10 — 15/i^^Mpc) l lDaddi et al.l 
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l2003l ; lGrazian et alj2006l) , comparable to the most luminous galax- 
ies in the local universe. 

In this paper we present a study of the first large sample of 
DRGs selected at bright infrared magnitudes (Kab < 21) in a 
contiguous area. We analyse the number counts and clustering and 
draw conclusions about their likely origin. Throughout this paper, 
we assume flm = 0.3, r2A = 0.7 and h = Ho/70 km s^^ Mpc^^. 



2 UKIDSS UDS EARLY DATA RELEASE 

2.1 Survey and Early Data Release 

The U KIRT Infrared Deep Sky Survey (UKIDSS - iLawrence eTal] 
is began observations in Spring 2005, using the Wide-Field 
Camera (WFCAM - Casali et al., in prep.) at the 3.8-m United- 
Kingdom InfraRed Telescope (UKIRT). Comprising 5 sub-surveys, 
it will take 7 years to complete and will cover a range of areas and 
depths. The deepest of these 5 sub-surveys, the Ultra-Deep Sur- 
vey (UDS) aims to cover 0.8 deg^ to a depth of Kab = 25.0, 
Hab = 25.4, Jab ~ 26.0 (5(t, point- source). It i s centred on the 
Suby u/XMM-Newton Deep Survey field (SXDS - ISekiguchi et all 
l2005h at 02''18™00^ -05°00'00" (J2000). 

Since February 10 2006, the UKIDSS Early Data Release 
(EDR) has been available to the E SQ com munitQ A full descrip- 
tion of this data release is given in lDve et al. (.2006.) . 

2.2 Image stacking and mosaics 

The stacking of the UDS EDR data was performed by our team 
using a slightly different recipe to the standard UKIDSS pipeline. 
Given the (relatively) small size of the field, it is possible to create 
a full mosaic before extracting catalogues, rather than merge cata- 
logues extracted from individual chips. Given the various jitter and 
offset sequences, this helps to optimise depth in overlap regions and 
to produce a more homogeneous final image. 

Each observation block consists of a 4-pt mosaic to tile the 
0.8 deg^ field, prod ucing 16 images each of ~ 15 minutes expo- 
sure (see lDveetal]r2006) . Individual reduced frames for each ob- 
servation block were extracted from the UKIDSS pipeline as the 
starting point for our final mosaiced stack (including astrometric 
and photometric solutions). We used the variance map produced 
by the pipeline to weight each frame before stacking and rescaled 
the pixel flux for each individual image using the pipeline zero- 
points. The stacking was carried out in a two-step pr ocess by us- 
ing th e SWarp software, an image resampling tool teertin et al.l 
l2002h . The final mosaiced images in the J- and if-band have the 
same pixel scale of 0.1342", with identical field centres and im- 
age scales to simplify catalogue extraction. The resulting images 
were visually inspected, and bad regions were masked out (areas 
around saturated stars, cosmetic problem areas, and low signal-to- 
noise borders). After masking, the usable area of this frame with 
uniform coverage is 0.62 deg^. We note that this is smaller than the 
expected 0.8 deg^ mosaic since the UDS field centre was moved 
by ~ 8 arcmin shortly after observations began (to allow the use of 
brighter guide stars). 

The image seeing measured from the mosaiced images is 
0.69" FWHM in K and 0.80" in J. The RMS accuracy of our 
astrometry is ~ 0.05" (i.e. < 1 pixel) and for our photometry 
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the R MS accuracy is < 2% in both J- and K- bands jDve et akl 
l2006b . From direct measurements of noise in a 2" aperture on the 
image we estimate 5a limiting magnitudes of Kab = 22.5 and 
JAi3=22.5. 



2.3 Catalogue extraction 

We found that the standard UKIDSS source detection software 
did not produce optimal catalogues for the UDS. We therefore 
produced a much improv ed catalogue for the E DR by using the 
SExtractor software (Bertin & AmoutsI 19961) . The if -band im- 
age was used as the source detection image, since this is measurably 
deeper for most galaxy colours. All Kab magnitudes quoted below 
are total magnitudes extracted using the SExtractor parameter 
MAG_AUTO, while all colour measurements are obtained from fixed 
2" aperture magnitudes. 

To optimise our catalogue extraction we performed a series 
of simulations to fine tune the SExtractor parameters. Artifi- 
cial point-like sources were added to the real Tf-band image using 
the observed PSF with FWHM =0.69" (rejecting regions contain- 
ing bright sources), and distributed with magnitudes in the range 
14 < Kab < 24. From the resulting new image a catalogue was ex- 
tracted using SExtractor and compared with the list of artificial 
source positions. This process was repeated 1000 times, and the re- 
sulting statistics allow us to estimate the catalogue completeness 
and the evolution of photometric errors. Using these simulations, 
SExtractor detection parameters were tuned to maximise com- 
pleteness at the noise-determined 5a depth of Kab = 22.5, while 
simultaneously minimising the number of spurious sources. While 
formally optimised for point-like sources, we note that these were 
close to optimal when we generated artificial sources using a sub- 
stantially more extended PSF (FWHM= 1.2"). 

Using these parameters, we extracted 78709 sources over 
0.62deg^ from the image, of which 34098 were determined to be 
unsaturated, unmasked and from regions of uniform coverage to 
Kab < 22.5. These form the basis of the analysis outlined below. 
Assuming the background noise is symmetric about zero, we can 
estimate the spurious fraction by extracting sources from the in- 
verted image and comparing with the number of sources extracted 
from the normal image. At our magnitude limit of Kab = 22.5 the 
fraction of spurious detections is found to be less than 1%, while 
the completeness level is above 70% (for point sources). 

3 SELECTION AND NUMBER COUNTS 
3.1 Selection of DRGs 

From the catalogue described above we selected objects using the 
(J — K)ab > 1.3 criteria. A visual inspection of each source was 
then required to remove spurious detections, which at these ex- 
treme colours was found to be a relatively large fraction 20%). 
The majority are caus ed by diffraction spikes and cross-talk im- 
ages toveetal.ll200^ and are easy to identify and reject. This 
leaves 369 DRGs at Kab < 21.2, which represents the largest sam- 
ple selected over a contiguous area. The surface density derived is 
n = 0.163 ± 0.009 arcmin"^. Figure[T]shows the (J — K) colour 
of these galaxies versus isT-band magnitude. The object shown by 
a star was classified as a point-like source in our global catalogue, 
and is confirmed to be a star after visual inspection. 
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Figure 1. Lower panel: (J — K)ab against Kab for sources in tlie UDS 
EDR field. Small points represent the full population of K-band selected 
objects, while larger points are those selected with (J — K)ab > 1-3 (and 
visually confirmed). The colour criteria, the magnitude limit at Kab = 
21.2 and the magnitude completeness at Kab = 20.7 are highlighted, as is 
the crude boundary between galaxies and the stellar locus at (J — K)ab = 
0. Upper panel: Errors as derived by SExtractor on the (J — K)ab 
colours (for display purposes only l/5th of the points are shown). Mean 
values are also displayed, as are errors derived from simulations. 
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Figure 2. A'-band differential number counts for our sample of DRGs. 
The errorbai's plotted are computed from poissonian small number statistics 
i Gehrels 1986). For comparison, we have overplotted the number counts 
for DRGs derived from the DXS EDR survey, with errors representing the 
field-to-field variance (see section IT4l - these points are slightl y shifted for 
display purposes), and from the literature from fainter samples j Labbe et"al] 



l2003tlF6rster Schreiber et alj2004rReddv et al.l2"00^lGrazian et alj20O6l) 
EROs number counts from the UDS/DXS SV sample are shown as well, for 
comparison. The plot inset shows number counts over a larger magnitude 
range. 



3.2 Photometric errors and contamination of DRG sample 

Since our sample is based on (J — K) colour selection it is vital 
to carefully consider the effects of photometric errors. Since most 
galaxies show substantially bluer colours (Figure[T]l we can expect 
the number density of DRGs to be artificially boosted at fainter 
magnitudes, as errors push objects above the (J — K)ab > 1.3 se- 
lection boundary. As a lower limit to this contamination we could 
use the photometric errors derived from SExtractor, and these 
are shown as a function of magnitude in Figure [T] Our experience 
suggests that analytically-determined errors from SExtractor 
are likely to be underestimates, so we use the mean photometric 
errors obtained from the simulations described in section [23] 

We used our simulated errors to estimate the contamination by 
randomising the real galaxy catalogue using Monte-Carlo simula- 
tions. For each object in our full catalogue, we allow the (J — K) 
colour to vary assuming a Gaussian distribution with a standard 
deviation corresponding to the chosen photometric error. We then 
re-select our catalogue using the (J — K)ab > 1-3 criteria, and 
repeat this process 1000 times. This should provide an approxi- 
mate upper limit on contamination, since we are randomising the 
observed galaxy catalogue (which already suffers from the effects 
of photometric errors). 

Defining the contamination fraction from the number of ob- 
jects scattered into our selection boundary minus those which are 
scattered out, our simulated source errors yield contamination frac- 
tions of (46.0 ± 3.8)% at the limiting magnitude of Kab < 21.2, 
falling to (16.8 ± 3.6)% at Kab < 20.7 (the estimated complete- 
ness limit; see Section [33). We note that the typical error on the 
colour is A(J — K) ab ^ 0.1 at Kab = 20.7. As shown in 
IConselice et al.l l l2006h . a slight change of the (J — K) colour se- 
lection does not have a major affect on the redshift distribution. 



Based on these values we conclude that our number counts 
and clustering measurements are reasonably robust at Kab < 20.7, 
but will become increasingly unreliable at fainter magnitudes. We 
will therefore adopt a limit of Kab ~ 20.7 for further study, which 
produces a sample of 239 DRGs. 



3.3 Number counts of DRGs 

Figure |2] shows the Tf-band differential number counts of our sam- 
ple of DRGs. The number counts indicate that our sample is com- 
plete up to approximately Kab — 20.7, after which the counts are 
clearly dropping. This defines our estimated completeness limit. 
Our simulations suggest that the contamination due to photometric 
errors will be ~ 16% at Kab < 20.7. We conclude the dominant 
source of error in our number counts will be Poisson counting er- 
rors (plotted) and cosmic variance (discussed in section l?!4t . 

At bright magnitudes (e.g. Kab < 20) our UKIDSS data are 
entirely unique, and no studies exist in the literature for compari- 
son. At fainter magnitudes, ou r counts are in very g ood agreement 
with the DRG counts from the lGrazian et al.l ( l2006h sample. They 
are also consistent with the numb er counts from the AE GIS survey 
(Foucaud et al. in prep.; see also IConselice et al]|2006l) . Combin- 
ing literature data with the present work we can examine the global 
shape of the DRG number counts over a very wide dynamic range 
(18.5 < Kab < 25.0). This strongly suggests a break feature in the 
slope at Kab ~ 20.5 which is an effect alr eady seen in the global 
K-band number counts (e.g. lGardner et airi 993). 

We note that the projected density of DRGs is approximately 
10 times lower than EROs, and approximately 100 times lower than 
the global galaxy counts at a given magnitude. 
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Figure 3. Distribution of our sample of DRGs on tlie sky, providing a visual 
impression of the large scale structure. Small dots represent the full K-band 
galaxy sample, while larger symbols are the DRGs with K^g < 20.7. 



3.4 Cosmic variance 

As a simple test of cosmic variance, and to investigate whether 
the UDS is an unusual field, we used the data available from the 
UKIDSS Deep Extragalactic Survey (DXS) to perform a compar- 
ison study. The DXS is the other deep extragalactic component 
of UKIDSS i Lawrence et al.ll2006l) . consisting of 4 fields with a 
7-year goal of observing 35 deg'^ to depths of Kab ~ 22.7 and 
Jab = 23.2. We used data from 3 fields observed in both J- and 
Tf-bands in the UKIDSS EDR, covering ~ 2900 arcmin^ ~ 4500 
arcmin^ and ~ 2900 arcmin^ respectively. While exposure times 
are similar to the UDS EDR, the observing conditions are gener- 
ally poorer. Direct comparison is also complicated by the different 
source extraction methods used by the DXS. 

We applied the same selection method described in sec- 
tion |3T| except that we did not visually inspect the samples. This 
selects 1523 objects in total, of which we estimate approximately 
20% are likely to be spurious (see section 13. it , with a similar 
fraction likely to be artificially boosted due to photometric errors 
at faint magnitudes (section |3.2| l. Since these errors are smaller 
than the errors in the DXS counts, for simplicity we opt not to 
make these corrections in our comparison with the UDS. We de- 
rived a median surface density of n = 0.176 ± 0.075 arcmin^^ at 
Kab ~ 21.2, in very good agreement with the UDS value. 

The resulting median counts from the 3 DXS samples are 
overplotted in figure |2l with errors representing the field-to-field 
RMS variance. The agreement with UDS is very good. Although 
no corrections were applied, this crude comparison suggests that 
the density of DRGs is stable and broadly consistent between fields. 



4 THE CLUSTERING OF DRGS 
4.1 Angular clustering 

In Figure[3l we display the distribution of our sample of DRGs on 
the sky. Visually the DRGs appear strongly clustered. As a more 
quantitative measure, we evaluate the 2-point angular correlation 
function, cutting our sample at Kab = 20.7 as before. 



T 1 — I M M 1 1 1 r — \ — I I I 1 1 1 r 




0.01 0.1 
0(deg) 



Figure 4. 2-point angular correlation function determined for our sam- 
ple of DRGs (full circle) and for the field galaxies (open circles) with 
Kab < 20.7 and their fitted power-laws . The plot inset presents the 
minimisation to estimate the best-fitted values of Ai^ and S for our sample 
of DRGs, with contours showing Icr and 3cr confidence levels. For compari- 
son, measurements from the literature for different samples are ove rplotted, 
assummg a slope of 5 = 0.8 teong et alj2006l : lGrazian et alj2006l) . 



To measure the angular correlation fu nction cij (6') and e s timate 
the related poissonian errors we used the iLandv & Szaiavl ( Il993h 
estimators. Figure |4] shows the correlation function derived from 
our DRG sample. The best fit for the angular correlation is assumed 
to be a power-law as: 

uj{6)=A^{6-' -Cs) 

with the amplitude at 1 degree — 3.11^ 3 x 10^^, the slope 
5 = 1.00 ± 0.10, and the integral constraint due to the limited area 
of the survey Ca =3.71 (determined over the unmasked area). 

For comparison we also determine the angular correlation 
function for the full K-selected sample to a limiting magnitude of 
Kab < 20.7, as shown in figure|4] The derived amplitude at 1 de- 
gree is yl„ =8. ll^j X lO"'' for a slope of 5 = 0.96 ± 0.05. Fixing 
the slope to the standard literature value (S = 0.8) we derived an 
amplitude of = 1.8ln i" x 10~ ^ which is in good agreement 
with the work o f lKong et al.N2006l) to the same magnitude limit. 

As sho wn in figure [4] by comparing with measurements from 
iKong et a P ( 2006) at the same magnitude limit of Kab < 20.7, the 
amplitude of the DRG sample is ~ 5 times higher than that of field 
galaxies but a factor of ~ 2 lower than EROs and BzK-selected 
galaxies. Interestingly, our derived clustering amplit ude is notably 
higher than measurements for DRGs obtained by iGrazian et al.l 
(2006) . This issue is explored further below. 

4.2 Spatial correlation lengths and biasing 

In order to compare the clustering of galaxy populations at different 
redshifts we must derive spatial correlation measurements. Assum- 
ing a redshift distribution for our sample we can derive the correla- 
tion length and a linear bias estimation using the relativistic Limber 
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equation dMagliocchetti et al.l200(]|) . The difficulty here is to have a 
realistic estimation of the redshift distribution of our sample. While 
Ipranx et aP l l2003h have designed the ( J — K )a b > 1.3 criteria 
to sele ct z > 2 galaxies. iGrazian et al.l ( |2006[) and lConselice et all 
( l2006h have shown that the redshift distribution is broad, and the 
fraction of z < 2 galaxies increases at bright magnitudes. 

As a preliminary investigation we assume a Gaussian form 
for the redshift distribution, fixing the mean and standard devia- 
tion to match recent studies of fainter DRGs. Analytical and fitted 
redshift distributions make no significant difference to our derived 
correlation lengths. The correlation length can then be evaluated if 
we fix the slope of the correlation fu nction to 7 = 1 + 5 = 2.0. 
In the DRG sample of iGrazian et alj ((2006), a mean redshift of 
z = 1.5 was found at a limit of Kab < 22.0. Using this mean 
redshift with ct = 0.5 (broadly matching their distribution) we de- 
termine rp = M-ltt o ft~^Mpc and 6 = 4.8ti J for o"'' DRGs. 
IConselice et al. suggest that very bright DRGs lie at even 

lower redshifts. If we adopt the observed redshift distribution of 
their spectroscopic sample we derive ro — ll.Otj 3 h'^Mpc and 
b = 4.0lj g, when the approximated gaussian form gives (with 
2 = 1.0 a nd g = 0.25) ro^ ll.ltyj h'^Mpc and 6 = 4.0t J:^. 

Grazian et alj {200^ also split their sample according to red- 
shift, and found ro = 7.4^4 9 h'^Mpc for DRGs at 1 < z < 2 with 
Kab < 22. Our sample is substantially brighter, so, even if the two 
samples are probing slightly different redshifts, our larger correla- 
tion length could be interpreted as evidence for luminosity segre- 
gation, e.g. we may b e sampling more l uminous galax ies which are 
more strongly biased jPollo et alj|2006h . We note that lOuadri et al.l 
( I2OO6I) find no relationship between K-magnitude and correlation 
length for K-selected galaxies, but instead find a strong relationship 
between colour and correlation length (with redder galaxies being 
more strongly clustered). This is not inconsistent with our findings, 
since we are studying segregation within a colour-selected sample. 

As shown in figure |4l our measurement is comparable to the 
clustering of EROs or BzK selected galaxies at the same limit- 
ing magnitude, which are likely to sample a similar redshift distri- 
bution to our DRGs, and hence show similar correlation lengths. 
lOverzier et alj i2003h compared the clustering of EROs with pre- 
dictions of Dark Matter halo models. Adopting the same technique, 
we can infer that bright DRGs are likely to be hosted by very mas- 
sive dark matter halos (possibly M > 10^'^ Mq), consistent with the 
potential progenitors of present-day L>L, elliptical galaxies. 

Finally we note that while the slope of our angular correlation 
function is steep, the fiducial value of 5 = 0.8 is only excluded by 
our sample at a la-level. Further data will be required to investigate 
whether this may be caused by a possibl e transition betwe en inner 
and outer halo occupation numbers (e.g. iOuchi et al.ll2005h . 



5 SUMMARY 

We have extracted a large sample of bright Distant Red Galaxies 
from the UKIDSS UDS EDR. Our catalogue contains 369 DRGs 
to a limiting magnitude of Kab =21.2, extracted over an area of 
0.62 deg^. The fainter Kab > 20.0 number counts are in good 
agreement with previous estimates, while at brighter magnitudes 
the sample is unique. Using simulations we determined that con- 
tamination due to photometric errors is below ~ 16% at an approx- 
imate completeness limit of Kab < 20.7. 

From this sample we extracted a sub-sample of 239 bright 
DRGs to a limit of Kab = 20.7. These bright DRGs appear 
highly clustered, and we determine a correlation length of ro — 



11 /i~^Mpc and a bias measurement b ~ 4.5, assuming the sam- 
ple lies at a mean redshift of z = 1.0 with a stand ard deviation of 
(7 = 0.25 (consistent with studies at similar depths - IConselice et al.l 
[2OO6). They appear more clustered than fainter samples of DRGs 
derived at these redshifts, which may be evidence for luminosity 
segregation, in agreement with biased galaxy formation scenarios. 
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